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Abstract 

We calculate the CP averaged branching ratios and CP-violating asymmetries for B® — ► gn, gn' 
and g'g' decays in the perturbative QCD (pQCD) approach here. The pQCD predictions for 
the CP-averaged branching ratios are Br(B® — ► gg) = (14.2+^°) x 10 -6 , Br(B® — > nn') = 
(I2.4t7 8 2 ) x 10~ 6 , and Br(B° s -> rfrf) = (9.2^ 5 9 3 ) x 10~ 6 , which agree well with those obtained 
by employing the QCD factorization approach and also be consistent with available experimental 
upper limits. The gluonic contributions are small in size: less than 7% for B s — ► gg and gg' 
decays, and around 18% for B s — ► n'g' decay. The CP-violating asymmetries for three decays 
are very small: less than 3% in magnitude. 
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Among various B — > MiM 2 decay channels ( here M, refers to the light pseudo-scalar 
or vector mesons ), the decays involving the isosinglet rj or rj' mesons in the final state are 
phenomenologically very interesting and have been studied extensively during the past 
decade because of the so-called Krj' puzzle or other special features [2, 0, H| • 

Motivated by the large number of B s production and decay events expected at the 
forthcoming LHC experiments, the studies about the B s meson decays become more 
attractive than ever before. Very recently, some two-body B s — > Miij^ decays, such as 
B s — > (tt, p, ou, <f>)rj^ decays have been studied in Refs. [!, @] in the perturbative QCD 
(pQCD ) factorization approach @, @, @]- In this paper, we would like to calculate the 
branching ratios and CP asymmetries for the three B° s — > rjr), rjr)' and rj'rj' decays by 



employing the low energy effective Hamiltonian [10J and the pQCD approach. Besides 
the usual factorizable contributions, we here are able to evaluate the non-factorizable and 
the annihilation contributions to these decays. 

On the experimental side, only the poor upper limit on Br(B° s — > r/rj) is available now 
[HI (upper limits at 90% C.L.): 

Br(B° s -> r]r]) < 1.5 x 1(T 3 , (1) 

Of course, this situation will be improved rapidly when LHC experiment starts to run at 
the end of 2007. 

This paper is organized as follows. In Sec. [H we calculate analytically the related 
Feynman diagrams and present the various decay amplitudes for the studied decay modes. 
In Sec. HIl we show the numerical results for the branching ratios and CP asymmetries 
of B® — > rf^rf'} decays. A short summary and some discussions are also included in this 
section. 



I. PERTURBATIVE CALCULATIONS 

Since the b quark is rather heavy we consider the B s meson at rest for simplicity. It 
is convenient to use light-cone coordinate {p + ,P~,Pt) to describe the meson's momenta: 
p± = (p° ±p 3 )/\/2 and pr = (p 1 ,^ 2 )- Using the light-cone coordinates the B s meson and 
the two final state meson momenta can be written as 

Mr Mr Mr 

A = ^|(1,1,0 T ), P 2 = -^(1,0,0 T ), P 3 = ^|(0,1,0 T ), (2) 

respectively, here the light meson masses have been neglected. Putting the light (anti-) 
quark momenta in B s , 7/ and rj mesons as ki, k 2 , and k 3 , respectively, we can choose 

k 1 = (x 1 Pf,0,k 1T ), k 2 = (x 2 P 2 + ,0,k 2T ), k 3 = (0,x 3 P 3 ~,k 3 T). (3) 

Then, after the integration over fcj", fcj, and k 3 , the decay amplitude for B s — > 777/ decay, 
for example, can be conceptually written as 

A(B S — » r]7]') ~ Jdxidx2dx 3 bidbib2db2b 3 db 3 

•Tr [C(t)<S> Bs (x 1 ,b 1 )<S> n ix2 1 b2)%(x 3 ,b 3 )H(x l ,b l ,t)S t (x t )e- s ^] , (4) 
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where k{ are momenta of light quarks included in each meson, term Tr denotes the trace 
over Dirac and color indices, C(t) is the Wilson coefficient evaluated at scale t, the function 
H(ki, k 2 , k 3 ,t) is the hard part and can be calculated perturbatively, the function $ M is 



the wave function, the function St(xi) describes the threshold resummation [12j which 
smears the end-point singularities on Xi, and the last term, e~ s ^\ is the Sudakov form 
factor which suppresses the soft dynamics effectively. We will calculate analytically the 
function H(xi,bi,t) for the considered decays in the first order in a s expansion and give 
the convoluted amplitudes in next section. 

For the two-body charmless B s meson decays, the related weak effective Hamiltonian 
H e ff can be written as fTT 



n 



eff 



G p 

7! 
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v uh v* q (cmo u m + cmo u M) - v tb v* £ c^) o^) 



i=3 



(5) 



where Cj(/i) are Wilson coefficients at the renormalization scale u and Oi are the four- 
fermion operators for the case of b — > q (q = d, s) transition T(| . For the Wilson 
coefficients Cj(//) (i — 1, . . . , 10), we will use the leading order (LO) expressions, although 



the next-to-leading order (NLO) results already exist in the literature lOf j . This is the 
consistent way to cancel the explicit fi dependence in the theoretical formulae. For the 
renormalization group evolution of the Wilson coefficients from higher scale to lower scale, 
we use the formulae as given in Ref.[l3[ directly. 



A. Decay amplitudes 

We firstly take B s — > r\r{ decay mode as an example, and then extend our study to 
B s — > rji] and r(vl decays. Similar to the — > 7i°i]^ decays in [BJ, there are 8 type 
diagrams contributing to the B s — > r/r] decays, as illustrated in Fig.l. We first calculate 
the usual factorizable diagrams (a) and (b). Operators 0^2,3,4,9,10 are (V — A)(V — A) 
currents, the sum of their amplitudes is given as 

rl poo 

F erj = 8TcCF'm I3s / dxidxs / bidbib 3 db 3 (j)B 3 (xi,bi) 
s Jo Jo 

X { [(1 + x 3 )rf(x 3 , 63) + (1 - 2x3)^(0^x3, 63) + ^(x 3 , b 3 ))] 

■a s {t\) h e (xi, x 3 , 61, 63) exp[-S ab (tl)} 
+2r v <f)%(x 3 , b 3 )a s {t 2 e )h e (x 3) x x , b 3 , h) exp[-S ab (tl)}} . (6) 

where r v = m^/rriB', Cp = 4/3 is a color factor. The explicit expressions of the function 
h e , the scales f e and the Sudakov factors S ab can be found Ref. [5|. The form factors of 
B s to 1] decay, F^ _> ' 7ss (0), can thus be extracted from the expression in Eq. 

The operators 5i6j7j8 have a structure of (V — A)(V + A). Some of these operators 
can contribute to the decay amplitude in a factorizable way, but others may contribute 
after making a Fierz transformation in order to get right flavor and color structure for 
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FIG. 1: Typical Feynman diagrams contributing to the B s — > ijrj' decays, where diagram (a) 
and (b) contribute to the B s —* rj form factor F { 



factorization to work. Such kinds of contributions can be written as 



7 pi 

eri 
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I dx\dx 3 / bidbib 3 db 3 (pB 3 (xi,bi 
o Jo 



(7) 



3 V 2m s m Bs 
x { [4>^(x 3 , b 3 ) + r v ((2 + x 3 )(f)P(x 3 , b 3 ) - x 3 (p^(x 3 , b 3 ))] 
■oi s {tl)h e (x u x 3 , b x , b 3 ) exp[-Sab(tl)] 
+ [xi0^(x 3 , b 3 ) - 2(xi - 1)^0^X3, 63)] 
■a s (t 2 e )h e (x 3 , x u 63, 61) exp[-S ab (t 2 e ))} . 



(8) 



For the non-factorizable diagrams 1(c) and 1(d), the corresponding decay amplitudes 
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can be written as 



M ev = — - — ■nC F m 4 Bs I dx 1 dx 2 dx 3 / b 1 db 1 b 2 db 2 4>b s ( x i, bi)4>%( x 2, h) 
-J Jo Jo 



x { [2x 3 r v (f)T(x 3 , bi) - x 3 (f) v (x 3 , bi)] 
■a s (t f )h f (x u x 2 , x 3 , b 2 ) exp[-S cd (t f )]} , 



M£ = 0, 

= -M ev 



For the non-factorizable annihilation diagrams 1(e) and 1(f), we find 

16a/6 f 1 f°° 

M ar] = — - — ixC F m A Ba I dx 1 dx 2 dx 3 / bidbib 2 db 2 </>B 3 (xi,b 1 ) 
<J B Jo Jo 

x { - {x 2 (p^(x 3 , b 2 )(j)^(x 2 , b 2 ) 

+r v r v > [(x 2 + x 3 + 2)0^ (x 2 , b 2 ) + (x 2 - x 3 )4\ ]l (x 2 , b 2 )] (p^(x 3 , b 2 ) 
+r v r v , [(x 2 - x 3 )(f)^(x 3 ,b 2 ) + (x 2 + x 3 - 2)<f)^(x 2 ,b 2 )] 4?(x 3 ,b 2 )} 

■a s (t 3 )h 3 (x 1 ,x 2 , x 3 , fei, b 2 ) exp[-S ef (t 3 )] 
+ {x 3 ^(x 3 ,b 2 )(l)^(x 2 ,b 2 ) 

+r v r v > [(x 2 + x 3 )(p^,(x 2 , b 2 ) + (x 3 - x 2 )4%,(x 2 , b 2 )] (p^(x 3 , b 2 ) 
+r r? ry [(x 3 - x 2 )(fr%(x 2 ,b 2 ) + (x 2 + x 3 )4%,(x 2 , b 2 )] (p^(x 3 ,b 2 )} 
•a s (t 4 f )h 4 f (x 1 ,x 2 , x 3 , 6i, b 2 ) exp[-S ef (t 4 f )]} , 

M S = -^—' R C F 'm A Bs j dx 1 dx 2 dx 3 \ b 1 db 1 b 2 db 2 (f)B s (x 1 ,b 1 ) 
<J Jo Jo 

x { [(x 3 - 2)r v ^,(x 2: b 2 )(^(x 3 , b 2 ) + (j)T(x 3 , b 2 )) - (x 2 - 2)r rf <j^{x 3 , b 2 ) 
((f)%(x 2 , b 2 ) + 4%,(x 2 , b 2 ))] ■ a s (t 3 )h 3 (xi,x 2 , x 3 , h, b 2 ) exp[-S ef (t 3 )] 
- [x 3 r v (f)$(x 2 , b 2 )((p^(x 3 , b 2 ) + 4%{x 3 , b 2 )) 
-x 2 r v ,(/)^(x 3 , fe 2 )(0j(x 2 , b 2 ) + 0j(x 2 , b 2 ))] 
■a s (tj)hj(xi,x 2 , x 3 , bi, b 2 ) exp[-S ef (tj)] } , 

Mar? — — o — ' K C F m Bs j dx\dx 2 dx 3 \ b 1 db 1 b 2 db 2 (p Bs (x 1 ,b 1 ) 

3 Jo Jo 

X {{ x 3^( x 3, b 2 )tf, (x 2 , b 2 ) 

+r r) r ri > [(x 2 + x 3 + 2)(p^(x 2 , b 2 ) + (x 3 - x 2 )<j^,(x 2 , b 2 )] <fr^(x 3 , b 2 ) 
+r 7? ry [(x 3 - x 2 )0j (x 3 , b 2 ) + (x 2 + x 3 - 2)^J (x 2 , b 2 )] <j%(x 3 , b 2 ) } 

■a s (t 3 f )h 3 f (x 1 ,x 2 , x 3 , 6i, b 2 ) exp[-S ef (t 3 )] 
- {x 2 (p^(x 3 , b 2 )(f)^(x 2 , b 2 ) 

+r v r v/ [(x 2 + x 3 )(pP(x 2 , b 2 ) + (x 2 - x 3 )(p^(x 2 , b 2 )] (p^(x 3 , b 2 ) 
+r ri r rt , [(x 2 - x 3 )(p^,(x 2 ,b 2 ) + (x 2 + x 3 )(f)^,(x 2 ,b 2 )] (p^(x 3 ,b 2 )} 
■a s (tj)h 4 f (x 1 ,x 2 , x 3 , 6i, b 2 ) exp[-S ef (t 4 f )]} . 



For the factorizable annihilation diagrams 1(g) and 1(h), we have 



/*1 /*oo 

F ar, = F P4 = -87iC F m Bs / dx 2 dx 3 / b 2 db 2 b 3 db 3 { [x 3 (j)^(x 3 , b 3 )(j)^{x 2 , b 2 ) 

Jo Jo 

+2r v r J ,,((x 3 + l)(p P (x 3 ,b 3 ) + (x 3 - l)(j> r (x 3 ,b 3 ))(j)^(x 2 ,b 2 )] 
■a s (t 3 e )h a (x 2 , x 3 , b 2 , b 3 ) exp[-S gh (t 3 e )] 

- [x 2 (p^(x 3 , b 3 )(f)$(x 2 , b 2 ) 

+2r v r v >((x 2 + l)(f>^(x 2 , b 2 ) + (x 2 - 1)0^, (x 2 , b 2 ))(p p (x 3 , b 3 )] 

■a s (t 4 e )h a (x 3 , x 2 , b 3 , b 2 ) exp[-S gh (t 4 e )]} (15) 

= -lQirC F m 4 Ba / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x { [x 3 r v ((j) P (x 3 , h) - $(x 3 , b 3 ))^{x 2 , b 2 ) + 2r n ,(f)f t {x 3 , b 3 )(j) P ,(x 2l b 2 )] 

■a s (tl)h a (x 2 , x 3 , b 2 , 63) exp[-S gh (tl)] 
+ [x 2 r n >{^{x 2: b 2 ) - 0j(x 2 , b 2 ))<fi^(x 3 , h) + 2r n ^,{x 2 , b 2 )<fi P (x 3 , b 3 )] 

■a s ( y t 4 )h a ( y x 3 ,x 2 ,b 3 ,b 2 )exp[-S gh (t 4 e )}} . (16) 

For the B s — > r)T)' decay, besides the Feynman diagrams as shown in Fig. 1 where 
the upper emitted meson is the 77', the Feynman diagrams obtained by exchanging the 
position of 77 and 77' also contribute to this decay mode. The corresponding expressions 
of amplitudes for new diagrams will be similar with those as given in Eqs. (loT THj) . since 
the 7] and 77' are all light pseudoscalar mesons and have the similar wave functions. The 
expressions of amplitudes for new diagrams can be obtained by the replacements 

</>$ < — * </>*/, </>% < — * 0j> 4^ < — ► r v < — * r v'- ( 17 ) 

For example, we find that: 

F F F , = -F F p } = -F P1 F p ? = F P2 (18) 

erf 1 er\i clij L ar\i L a r\' 1 arj ' 1 arf 1 at) " \ ±a ) 

Before we write down the complete decay amplitude for the studied decay modes, we 
firstly give a brief discussion about the 77 — 77' mixing and the gluonic component of the 
77' meson. There exist two popular mixing basis for 77 — 77' system, the octet-singlet and 
the quark flavor basis, in literature. Here we use the SU(3) F octet-singlet basis with the 
two mixing angle (01,0s) scheme []3] to describe the mixing of 77 and 77' mesons. In the 
numerical calculations, we will use the following mixing parameters 14j 

8 = -21.2°, 0x = -9.2°, h = 1.17U, f 8 = 1.26U (19) 

In this paper, we firstly take 77 and rj' as a linear combination of light quark pairs 
mm, dd and ss, and then estimate the possible gluonic contributions to B s — > 
decays by using the formulae as presented in Ref. [4|]. We found that the possible gluonic 
contributions are indeed small. 
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B. Complete decay amplitudes 



For — > r/r]' decay, by combining the contributions from different diagrams, the total 
decay amplitude can be written as 



M{rrf) = {F ev F 2 f^ + F ev ,F^) 



Zu[C 1 + -C 2 



— 6 ^C*3 + — C 5 — -C 6 + -C-j + -C 8 — -C 9 — -Cio 

- (F ev F 2 f°, + F er) ,F 2 f°) 

(75 2 f 1 f f 

3^3 + 3C4 - 2C 5 - -C7 6 - -CV - -C 8 + -C 9 - -C 10 

+ (M e , + M e ,0 F 2 F 2 ' £ U C 2 - 6 (c 3 + 2C 4 - ic 9 + ^C w 

- [M^F^ + M er) iF[F 2 ] & (c 4 - ^C 10 ) 

- « 2 + < 2 ') ^6 (2C 6 + ic 8 

- [M^FiF 2 + M%F[F 2 ) & (c 6 - ±<7 8 
+ (M ar) + M ar) ,) F X F[ L<7 2 - 6 f C 3 + 2C 4 - ic 9 + \c w 



- (M ar) + M w ) F 2 F& (c A - l -C^ 
~ (MS + M* ) F 2 F£ t (c 5 - l -C 7 

- « 2 + < 2 ' ) FiKit (2a + \c x 
~ « 2 + < 2 ' ) F 2 F^ t (c 6 - l -C, 
-f B . ■ K 2 + F%) F 2 F£ t Qc 5 + C, 



—C7 — — Cf. 

6 2 * 



(20) 
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where £ u = V* b V us and £ t = V t lV tS} and the relevant mixing parameters and decay con- 
stants are 

F x = y|cos0 8 - J|ffln0i, F 2 = -y|sin0 8 + y|cos0!, (21) 

F{ = Y^sin 8 + cos 1; F 2 ' = - J| sin 8 + cos 1; (22) 

# = Acos0 8 - A sin 1; /* = -^|cos0 8 - 4sin0i, (23) 

# = -| s in0 8 + A cos ^ /• =-^^^ + ^008^ (24) 

Similarly, the decay amplitudes for B® — > 7777 and 5° — > 7/7/ decay can be obtained 
easily from Eq.( j20l) by the following replacements 

& f'n^fi, ft\ Fi(0i, 08)^^(01,08); F 2 (0 1) 8 )^F 2 '(0 1 ,0 8 ). (25) 

Note that the contributions from the possible gluonic component of rf meson have not 
been included here. 



II. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we will calculate the CP-averaged branching ratios and CP violating 
asymmetries for those considered decay modes. The input parameters and the wave 
functions to be used are given in Appendix |A] In numerical calculations, central values 
of input parameters will be used implicitly unless otherwise stated. 

Using the decay amplitudes obtained in last section, it is straightforward to calculate 
the branching ratios. By employing the two mixing angle scheme of 77 — v( system and 
using the mixing parameters as given in Eq. (TI91) . one finds the CP-averaged branching 
ratios for the considered three decays as follows 

Br( B° s -> 7777) = [14.2+^)^^)] * 10 ~ 6 > ( 2 6) 
Br{ B° s -> W) = [l2-4l^K)™(mn] x 10" 6 , (27) 
Br{ B° s -> 77V) = [9.2t|°(o; 6 )^ 5 5 (mr)] x 10~ 6 , (28) 

where the main errors are induced by the uncertainties of Ub = 0.50 ± 0.05 GeV, and 
m* s = [1.49 — 2.38] GeV (corresponding to m s = 130 ±30 MeV), respectively. The above 
pQCD predictions agree well with those obtained in the QCD facterization approach [2J. 

As for the gluonic contributions, we follow the same procedure as being used in Ref. [4j 
to include the possible gluonic contributions to the B s — > 77 W transition form factors 
F [ and found that the gluonic contributions to the branching ratios are less than 
3% for B — > 7777 decay, ~ 7% for B — > rjr) 1 decay, and around 18% for B — > tj't] 1 decay. 
The central values of the pQCD predictions for B s — > rf''rf^ decays after the inclusion of 
possible gluonic contributions are the following 

Br(B° s -> 7777) = [lS.T^t(cj b )tff(rn^)] x lO" 6 , 
Br( B° s -> 7777') = [ll.6tiiK)^ 6 7 8 (mr)] x 10- 6 , 
Br{ B° s -> 77V) = [I0.8±i:l(a7 6 )±| 6 2 2 (mg-* )] x 10~ 6 . (29) 
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Now we turn to the evaluations of the CP-violating asymmetries of B s — > rf^rf^ decays 
in pQCD approach. For B° s meson decays, a non-zero ratio (AT/T)b s is expected in the 



SM 0, M ■ For B s — » 77M77M decays, three quantities to describe the CP violation can 



be defined as follows 1161: 



A dir _ \^cp\ 2 ~ 1 A mix _ 2/m(A C p) 2Re(X CP ) , , 

Acp ~ i + \x CP r CP ~ l + iw Ars ~ i + iA CP r (30) 



with 



V t tV ts (f\H eff \B°) ^ (f\H eff \B°) 
CP Vf V tb VMHe ff \B°) ~ (f\H eff \B°Y { > 

in a very good approximation. Here Aq P and A™™ means the direct and mixing-induced 
CP violation respectively, while the third term *4.Ar s is related to the presence of a non- 
negligible Ar s . By using the mixing parameters in Eq. ffl9|) and the input parameters as 
given in Appendix A, one found the pQCD predictions for Aqp, Aqp and Hf 

Af P (B° s - m) = [-0.2±0.1( 7 )±0.lM^(mr)]xlO- 2 , 

A&W^rrf) = [+0.6^(7) ± O.IK) ± 0.3(mn] xlO" 2 , 

A d MB° s - V'V') = [-°- 8 -ai(7) ± 0.1(0*) ± 0.7K')] x 10" 2 , (32) 

Atp (B° s -> 7777) = [-0.3 ± 0.1(7) ± 0.2(w 6 ) ± 0.5 W)] x 10" 2 , 

^op (-8° W) = h - 8 ± 0-2(7) ± 0.1K) ± 0.2(mg")j x 1(T 2 , 

A^(B° S -> r/r/) = [+1.8l°J( 7 ) ± 0.0(u; 6 )l°j(mr)] x 10" 2 , (33) 

^Ar s (w) « A A r s (w') « ^Ar. W) « 1, (34) 

where the dominant errors come from the variations of CKM angle 7 = 60° ± 20°, cj& = 
0.50 ± 0.05 GeV and m v sS = [1.49 - 2.38] GeV ( corresponding to m s = 130 ± 30 MeV), 
respectively. It is easy to see that both the direct and mixing-induced CP violations of the 
considered B s decays are very small in magnitude, and thus almost impossible to measure 
them even in the LHC experiments. The above pQCD predictions are also consistent with 
the QCDF predictions 0,3. 

In short, we calculated the branching ratios and CP-violating asymmetries of B® — > 7777, 
7777' and 7/7/ decays at the leading order by using the pQCD factorization approach. Besides 
the usual factorizable diagrams, the non-factorizable and annihilation diagrams are also 
calculated analytically in the pQCD approach. From our calculations and phenomeno- 
logical analysis, we found the following results: 

• Using the two mixing angle scheme, the pQCD predictions for the CP-averaged 
branching ratios are 

Br(B° s -> 7777) = (14.2^ 8 5 °) x 1(T 6 , 
Br(B° s -> 7777') = (12.4^ 8 2 ) x 1(T 6 , 

Br(B° s - r/V) = (9.2™) x 1(T 6 , (35) 

where the various errors as specified previously have been added in quadrature. The 
pQCD predictions for the three decay channels agree well with those obtained by 
employing the QCDF approach. 
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• The gluonic contributions are small in size: less than 7% for B s — > 7777 and 7777' 
decays, and around 18% for B s — > rj'r]' decay. 

• The direct and mixing-induced CP violations of the considered three decay modes 
are very small: less than 3% in magnitude. 

Note added: After completion of this paper, the paper in Ref.[l8| appeared, and where 
a systematic study for the B s — > M1M2 decays in the pQCD factorization approach has 
been done. Since different mixing-scheme of 77 — rf system have been used, the explicit 
expressions of the decay amplitudes of the relevant decays are different in these two papers, 
but the numerical predictions for branching ratios and CP violations agree well with each 
other. The possible gluonic contributions are estimated here. 
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APPENDIX A: INPUT PARAMETERS AND WAVE FUNCTIONS 

In this Appendix we show the input parameters and the light meson wave functions to 
be used in the numerical calculations. 

The masses, decay constants, QCD scale and B° s meson lifetime are 



(/=4) 



A1£T 

MS 



250MeV, U 



1.4GeV, 
5.37GeV, 



m s = 
M w 



130MeV, 
130MeV, 



fx 



-- 230MeV, 
160MeV, 



80.41GeV,r B o = 1.46 x 10~ 12 s 



(Al) 



For the CKM matrix elements, here we adopt the Wolfenstein parametrization for the 
CKM matrix, and take A = 0.2272, A = 0.818, p = 0.221 and 77 = 0.340 [ll|. 
For the B meson wave function, we adopt the model 



!>b s (x,5) = N Bs x 2 (l - x) 2 exp 



2ul 2 ' 



(A2) 



where Uf, is a free parameter and we take Ub = 0.50 ± 0.05 GeV in numerical calculations, 



and Nb s = 63.67 is the normalization factor for Ub = 0.50. 

For the distribution amplitudes (f>^ _, cj)^ and </>^ _ , we utilize the result from the 
light-cone sum rule 13] including twist-3 contribution. For the corresponding Gegenbauer 
moments and relevant input parameters, we here use al dd 



0.115, al" 



-0.015, p Vdd 



m w /mQ, 773 = 0.015 and = —3.0. We also assume that the wave function of uu is the 
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same as the wave function of dd [3| . For the wave function of the ss components, we also 
use the same form as dd but with m s Q s and f y instead of and f x , respectively: 



/« = /*, f y = \jyi-fi (A3) 

These values are translated to the values in the two mixing angle method: 

fx = 152.1MeV, / 8 = 163.8MeV, 

0i = -9.2°, fl 8 = -21.2°. (A4) 
The parameters mf, (i = r}^r m \, r] sS ) are defined as: 

vm^) = w = rnl 2M 2 K - ml 

m ° - m °-(m u + m d y m ° - (2m.) " (A5) 
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